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In a recent work (Phys.Rev.C84, 044321, 2011) M.J. Ermamatov and P.R. Fraser have studied 
rotational and vibrational excited states of axially symmetric nuclei within the Bohr Hamiltonian 
with different mass parameters. However, the energy formula that the authors have used contains 
some inaccuracies. So, the numerical results they obtained seem to be controversial. In this paper, we 
revisit all calculations related to this problem and determine the appropriate formula for the energy 
spectrum. Moreover, in order to improve such calculations, we reconsider this problem within the 
framework of the deformation dependent mass formalism. Also, unlike the work of Bonatsos et al. 
(Phys.Rev.C83, 044321, 2011) in which the mass parameter has not been considered, we will show 
the importance of this parameter and its effect on numerical predictions. 
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I. INTRODUCTION 

Thanks to its relatively simple structure the Bohr 
Hamiltonian [ij continues to play an undeniable role 
in the study of nuclear structure within collective mod¬ 
els in competition with more sophisticated meth¬ 

ods such as C^asiparticle Random Phase Approxima¬ 
tion (QRPA) @, g and Interacting Boson Model (IBM) 
g. Also, its advantage in respect to these microscopic 
methods resides in its ability to provide collective states 
eigenenergies and corresponding wave functions of nu¬ 
clei in analytical form. So far, the Bohr Hamiltonian 
has been widely used with a constant mass parameter 
gsi- Recently, this assumption has been reexamined 
in the framework of Deformation Dependent Mass For¬ 
malism (DDMF) [13) El] emphasizing the mass tensor of 
the collective Hamiltonian cannot be taken as a constant 
but it has to depend on the collective coordinates. Such 
a formalism allows enhancing the precision of numerical 
calculations of nuclear characteristics. Moreover, Jolos et 
al. [l6l4l9l| have shown this mass parameter should split 
into ground state band, /3-band and 7 -band coefficients 
for deformed nuclei. Each coefficient is set to its average 
value over the wave function of the corresponding band 
state. Following this later procedure, M.J. Ermamatov 
et al. have studied rotational and vibrational spectra 
of axially symmetric nuclei [ 2 ^. Their calculations have 
been based on an analytical energy formula that the au¬ 
thors claimed they have obtained in a previous work . 
However, the used formula in [^ together with the corre¬ 
sponding wave functions were inaccurate as we will show 
in this paper. Therefore, the calculated transition rates 
by the same authors are also questionable. Besides, the 
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Bohr Hamiltonian’s dependence on two separable collec¬ 
tive coordinates fd and 7 where P also represents nuclear 
shape deformation enables one to choose nuclear collec¬ 
tive potentials as a sum of two separate terms, namely: 
a /3-potential U(/3) and a 7 -term ^( 7 ). In the present 
paper where we revisit the M.J. Ermamatov et al. work 
[ 20 | with the purpose to improve their calculations within 
DDMF, the potential term V{ 0) is chosen to be equal 
to Davidson potential (2^ as in [ 2 ^ and the 7 -potential 
V ( 7 ) taken to be equal to the harmonic oscillator. Such 
a problem has been solved in but with equal mass 
coefficients by means of superwmmetric quantum me¬ 
chanical method (SUSYQM) [2^[23|- Furthermore, we 
will display the essential role played by the mass parame¬ 
ter in the evaluation of nuclear characteristics unlike the 
Bonatsos et al. work [iM3 in which this parameter has 
been hidden. Thus, the eigenenrgies formula and the cor¬ 
responding wave functions are derived by means of the 
asymptotic iteration method (AIM) [ 2 ^. This method 
has proved to be a useful tool when dealing with physical 
problems involving Schrbdinger type equations |26l - l28l |. 
This paper is organized as follows : In Section H the 
position-dependent mass formalism is briefly described. 
In section HI, we propose the Bohr Hamiltonian with 
three different mass coefficients, that we use in Section 
IV in accordance with deformation-dependent mass for¬ 
malism. The exact separation of the Bohr hamiltonian 
in the case of axially symmetric prolate deformed nu¬ 
clei and the solutions of angular equation are achieved 
in section V. The radial equation is given in Section VI. 
Analytical expressions for the energy levels and excited- 
state wave functions are presented in Sections VH and 
VHI respectively, while the B{E2) transition probabili¬ 
ties are given in the Section IX. Finally, the section X is 
devoted to the numerical calculations for energy spectra 
and B{E2) transition probabilities with their compar¬ 
isons with experimental data and the available IBM ones, 
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while Section XI contains the conclusion. An overview of 
the asymptotic iteration method is given in Appendix A. 
While in Appendix B, we give the used formulas for the 
calculations of B{E2). 


II. POSITION-DEPENDENT MASS 
FORMALISM 


The general form of the Hamiltonian with effective 
mass depending on position has been originally intro¬ 
duced by Von Roos 


H =-(x)Vm” 

4 1 


-I-(x)Vm” +V{x), (1) 


where V is the relevant potential and the parameters 
S', k', X' are constrained by the condition S'+k'+X' = — 1. 
Assuming a position dependent mass of the form [s^ 


m[x) = moM{x),M{x) = 7773 ^,/(a;) = l+g{x), 




( 2 ) 


where mg is a constant mass and M (x) is a dimensionless 
rosition-dependent mass, the Hamiltonian o becomes 


iJ = -—[/^(x)V/«V/^ 


4mo L 


+ /^(x)Vrvr +Vix), (3) 


with (5 -|-k-|-A = 2. It is known that this Hamiltonian 
can be put into the form 


^ = -^v7Rv/(x)Vv7R + K//(x), (4) 


with 


2 mo 


Veffix) = Vix) + £- 1(1 - d - X)f{x)vy{x) 


+ {l-S){l-X){S7f{x)r 


(5) 


Where is a collective variable and TT 2 p, is an operator 
of the conjugate momentum. In the intrinsic frame we 
obtain from Eq. 


H = - 


1 d 04 d 


4H(/3,7) W'* 9/3 ^ /32 sin 37 97 


^ sin 37 ^ 


dy 


_ Qk 

kit2,3 sin'(7- 


1 9 q4 9 
ap 


1 


sin 37 dj 


fc=l,2,3 


Q 


dl 4/32 ^ sin 2 (^ -Ink) j 

(7) 


4i?(/3,7) 


For small amplitudes of 7 -vibration around 7 = 0 and /3- 
vibration around /3 = /3o 0, the collective coordinates 

could be considered as separable in the axial symmetry 
nuclei case. Thus, we can consider three separable states 
of nuclei, namely : the ground state, the j3 and 7 vibra¬ 
tional states. Each one of these states will have its own 
mass parameter equal to its average value over the wave 
function of the state under consideration : 


1. The ground state mass parameter 

{g.s.\ g.s.) = Brot ( 8 ) 

where we consider the ground state rotational 
band; 

2 . the 7 -mass parameter 

(7 I H(/?,7) I 7 ) = H 7 (9) 

where we consider 7 -vibrational state; 

3. The /3-mass parameter 

(/3|i3(/3,7) |/3) = i3/3 (10) 

where we consider /3-vibrational state. 

The procedure described above assumes the use of pro¬ 
jection operators. Using Eqs. (IMTOl) . we obtain from Eq. 
([T]) the following Hamiltonian 


where S and A are free parameters. 


III. BOHR HAMILTONIAN WITH MASS 
COEFFICIENTS 

In the laboratory frame, the Bohr Hamiltonian can be 

written as [13 

( 6 ) 


H = - 




_I 1 1 1 9 „■ o„, 9 

2(* I 5 I i) I 9/3 + 9 ^^ 1^3797 


1 

w 


Q 


V - 




+m7) (11) 


where i = g.s., /3 or 7 band depending on which state 
is considered. In the case of a small axially symmetric 
deformation of nuclei, the Bohr Hamiltonian with three 
different mass coefficients can be written as [13 


H — Hrot + Hj + Hp 


( 12 ) 
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where 


in the form [s^ 


Hrnt — 




&BrotP' 




(13) 


1 d d Ql Vjj) 

25^/32 jd-f^d-f 2B^ 4 / 32^2 + ^2 


(14) 


and 

Hp 


( I 2 15 2 1 5 \ 

2 \Bp ^ p dp ^ B0 p dp) ^ 

(15) 


IV. BOHR HAMILTONIAN WITH DIFFERENT 
DEFORMATION-DEPENDENT MASS 
PARAMETERS 


fci.3 sin^T- i^fc) 

^ 7;iQi + QI + QI) + Qs ( -ql (20) 

3 \sin 7 3/ 

In the same context, we consider a wave function of the 
form 

4/(/3, 7 , 0,) = Fr,^L{P)Vu,,Kh)'Dij^M (21) 

where 25(0^) are Wigner functions of the Euler angles 
9i(i = 1,2,3), and L is the total angular momentum, 
where M and K are the eigenvalues of the projections 
of angular momentum on the laboratory-fixed z-axis and 
the body-fixed z'-axis, respectively. As a result, Eq. m 
can be approximately separated into three equations: 


To construct a Bohr Hamiltonian with a mass depend¬ 
ing on the deformation coordinate /3, in accordance with 
the formalism described in section II, 


I -So I i) 


(16) 


we have to follow the procedure in [I 3 . Since the defor¬ 
mation function / depends only on the radial coordinate 
/3, only the /3 part of the resulting equation will be af¬ 
fected. 

The final result reads [3 


2 (i I So M) 

P 


_ ( _ R'i f rj 

j\ V J J 


+ 


E 

fc=l,2,3 ' 


Ql 


Itt/c) 




_ F _a_ • 

sin 37 d'y ' d") 

+ = (17) 


with. 


P ( Pf 9 pin 9 /7 , Z!a 
2{i I So 7) V P'^dp^^dp'^^ /32 

+ i(l - 5 - A)/ VV + (i - ^)(i - A)(v/)^) 


+ V{P) 


S„,L(/3)=SS„^i(/3) 


( 22 ) 


tP f 1 d d 1 \ 

2Sq y sin 37 dj dj 4 sin^ 7 j 


+ W{j) 


Vn^,Kh) = 


(23) 


and 

^((^2 _ g 2 ^p^_^( 0 ,) = (24) 


+ (^-<5)(^-A)(v/)^) (18) 


The eigenvalues of the rotational part equation (|24)) are 
easily obtained since is the quadratic casimir operator 

of 0(3) and Qi is the projection of the angular momen¬ 
tum on the z-axis, 


V. SEPARATION OF THE BOHR 
HAMILTONIAN OR AXIALLY SYMMETRIC 
PROLATE DEFORMED NUCLEI 

Exact separation of the variables /3 and 'y may be 
achieved when the potential is chosen as in [l^ : 

1/(/3,7) = S(/3) + EiE(7) (19) 

where the potential 1 ^( 7 ) has a minimum around 7 = 0 . 
Then, one can write the angular momentum of Eq. m 


Note that Eq. (l23l) for 7 r; 0 can be treated as in 

Ref. M- 

For the 7 -part, we use a harmonic oscillator potential 

M 

11^(7) = \(PtC^)P (26) 

where /3o denotes the position of the minimum of the 
potential in /3 and is a free parameter. In this case. 
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Eq. (1231) transforms into the usual harmonic oscillator 
equation 


ti^ (l d d l\ 


+ 2 (^ 0 ^ 7 ) 7 ^ 


Vn^,K{'l) = AlJn^^Kh), (27) 


To solve this equation through AIM, we propose the fol¬ 
lowing ansatz for the 7 -part eigenvectors rjn^ 

»7n^,K(7) = (28) 

with g = ^ For this form of the angular wave 

^0 By Cy 

function, the 7 -part equation (HZl) reduces to a standard 
form given in the Appendix Eq. inj. According to the 
AIM procedure, the eigenvalues are calculated by means 
of the termination condition Eq. (IA5I) and the recurrence 
relations Eq. (TAdl) . hence one can derive the generalized 
form of the eigenvalues, 

A = -^2n^-I-= 0,1,2,..., (29) 

By inserting hj = in Eq. where is the 

quantum number related to 7 -oscillations, one obtains 

+ = 0,1,2,..., (30) 

g 

As a result, we found, 



(lEl 

\g 


{uj -(-1) -f 


1 Bp 
3 Brot 


(l{L + 1) 



(31) 


The allowed bands are characterized by 


n-y = 0, K = 0] 

Tlry =1, K = ±2; 

717 = 2, A: = 0,±4; ... (32) 


In the standard case of constant mass where B^ = Bp = 
Brot = 1 and /i = 1, our formula Eq. (1311) matches up 
with Eq. (41) of Ref [ij. In [T^, the coefficient of 7 ^ 
in u(j) is equal to (3c)^, compared to (/SqCj) Eq. (1^ 
used in this work. 

The eigenfunctions corresponding to eigenvalues (l30l) 
are obtained in terms of confluent hypergeometric func¬ 
tion, 


axially symmetric prolate deformed nuclei can be written 
as: 


r]n^,K 


= Nr, 




(?) 


K /2 

where represents the Laguerre polynomial and 

the normalization constant, determined from the normal¬ 
ization condition 



’7n.„K(7)|sin37|d7= 1 


(35) 


In the case of small 7 vibration, we can write | sin 87 ! ~ 
1 37 1, then the integral Eq. (I55|) is easily calculated by 
using Eq. (8.980) of [s^. This leads to 


IV, 


ny,K — 


? -l-|i<'/ 2 |__ 

3^ Tifi^ -4 |A:/2| -4 1) 


(36) 


The normalization constants for the {n^, K) = (0,0) and 
{rij, K) = (1, 2) states are found to be ^ iVf 2 = 

respectively, then = g. This result will be used 

to calculate the B{E2) values in 7 —ground and 7 —>■ 
P transitions (AAT = 2 ). 


VI. THE RADIAL SCHRODINGER EQUATION 


The /3-vibrations states of deformed nuclei with mass 
parameter are determined by the solution of the radial 
Schrodinger equation 


2 




7 ' 
/\2 


\Bp B^J\ 8 4 /3 / 


F 


-^AE + AE-K//E = 0 (37) 


with 

1 /I 


-{i-s-x)fr 


-y.-yyyf'r- 


(38) 


Setting a standard transformation of the radial wave 
function 


Fn,L(/3)=/3-(l+^^/^-)i?„,L(/3) 


(39) 


rn„K(7) = iEi( - 1 + y) (33) 

where is a normalization constant. According to 

the relation between hypergeometric functions and the 
Laguerre polynomials, the 7 angular wave functions for 


we get 




R 


+ 2UeffR=‘^ER (40) 
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where 


TT - I ^ /i I 

Ueff- ;j 2^-//+2 R (l+sJ 


Bp\f + /3ff , B0 f 


2/3" 


(41) 


In the frame without mass parameters, we can reduce the 
first three terms in the Eq. m by {vt^vd'r- 


VII. THE EFFECTIVE POTENTIAL AND 
ENERGY LEVELS 


As in Ermamatov et al. work 
calculations the Davidson potential 



we use in our 


VW) 



(42) 


where Vq represent the depth of the minimum, located at 

/3o. 

According to the specific form of the potential Eq. 
(HU), we are also going to consider for the deformation 
function the special form 

/(/3) = l + a/3^, a << 1 (43) 


Eor this form of the radial wave function, the Eq. HQl) 
reads. 


Xnpiy) = - 


—a 


1 + 4/) + ay{3 + 2-^ + 4p -h 4i/) 


2y{l + ay) 


x'npiy) 


2{p + v){\ + 2 ^ + 2 z/ + 2 p) + 1 + - I 

42/(1 + ay) 


k2 


Xnp{y) 

(48) 


The first and the second terms in the right hand side of 
Eq. (H5)) represent Aq and sq of Eq. (ED respectively. 
After calculating A„ and s„, by means of the recurrence 
relations of Eq. dH), we get the generalized formula of 
the radial energy spectrum from the roots of the termi¬ 
nation condition of Eq. (lASp 

„ /, a , Bg . 

Enpn^LK = [ko + —{2-\- — — \-2p + 2q + pq) 

-f 2u(2 p -\- 4o7r^^ (^9) 

where np is the principal quantum number of /3 vibra¬ 
tions, and 


q = qn^ (A, K) = ^/l + 4A:_2 


Inserting these forms for the potential and the deforma¬ 
tion function in Eq. m one gets 


with 


i, ® 

k2 = — 


2Ueff = k2P^ + ko + ^ 


(' + f)(4Dl-2^)(l-2A) 


(44) 


-b (l-(5-A)) -b 2 -^A 


,5/3 


i, ® 

fco =2 




_4^ 

«2 


(45) 


where gp = 

To solve the radial equation Eq. (TOl) through the 
asymptotic iteration method (AIM) [^, one needs the 
following parametrization 

Rnphiy) = 2/^(1 + ay)''XnpL{y), 2 / = /3^ (46) 

where 

1 


p=-[l + ^l + Ak- 2 ) 


2 




fco 2Bft 


1 / 2 - 


ah? 


E 


(47) 


p = p„jL,K) = + (50) 

The quantities ^ 2 , fco, ^-2 are given by Eq. (145L 
where A is the eigenvalue of the 7 —vibrational part of 
the Hamiltonian for axially symmetric prolate deformed 
nuclei. In the numerical results part of the paper, the 
energies are normalized to the first excited state. So, the 
results depend on six parameters Bp/B^, Bj/B^ot, 9, 9 / 3 , 
a and /3o. 

A few interesting low-lying bands are classified by the 
quantum numbers np, and A", such as the ground 
state band (g.s.) with np = 0, n.y = 0, K = 0, the 
/3—band with np = 1, = 0, K = 0, and the 7 —band 

with np = 0, = 1, K = 2. 


A. Special case 1: Without mass coefficients 


If we assume Bp = By = B^ot = 1, one gets from 

Eq.dlSl) 


= a^ 


kn = a 


(I - (5 - A) -b (I - 25) (I - 2A) -b 6 -b A 




(I - 5 - A) -b 8 -b 2A 


-4Eo 


fc _2 = 2 -b A -b 2Vb/3o (51) 

Thus, the energy spectrum formula Eq. (1491) is identical 
to Eq. (82) of Ref. 0 obtained by means of su persy ni- 
metric quantum mechanical method (SUSYQM) [2j,[^|. 
The slight difference between our coefficients fc 2 , and 
k _2 and those of Ref. [l^ comes from the adopted ex¬ 
pression of Davidson potential. 
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B. Special case 2: No dependence of the mass on VIII. EXCITED-STATE WAVE FUNCTIONS 

the deformation 


If a = 0, the dependence of the mass on the deforma¬ 
tion is canceled, then one obtains from Eq. (HSl) 



In this case, the energy spectrum becomes 
En/in^LK = 2 ^ + 2^/3 + -qn^{L,K))'j 

^ (53) 

For axially symmetric prolate deformed nuclei, the en¬ 
ergy formula reads 

/ y2 , 2. _\ 

Enpn^LK = W ~ 

(54) 

with 

(L, A) = (l + 1^) + + 2gp (55) 

and 

Note that Eq. (l5^ represents the correct formula of the 
energy spectrum, compared to Eq.(ll) given in Ref. [^ . 
where the mass parameter term is missed in the analogue 
formula of Eq. ([55]) . 

It is also worth to note that, in this case, Eq. (03 
reduces to standard confluent hypergeometric equation 
which can be converted to Laguerre differential equation. 
The resolution of such a problem is carried out in Section 
Vlll.b. 


C. Special case 3: Standard case 


For 7 -unstable nuclei, in the limit case of a = 0 and 
Bp = Bj = Brot, our formula Eq. reduces to 

+ 2n;3 + ^J^+A + 2gp^ - 2 Vo (57) 

with 


A = t{t + 3) (58) 

and T = L/2 is the seniority quantum number. This 
formula is similar to the energy spectrum Eq. (80) in 
Ref. m. 


The used wave functions in our calculations are given 

by 

'^{13,^,Oi) = Rnf),LiP)Vn^,Kh)E>MKi^i) ( 59 ) 

The radial function Rnft,L{P) corresponds to the 
eigenstate of Eq. (SOD, 9n^., K ( 7 ) is given by Eq. (IMl) and 
the symmetries eigenfunctions of the angular momentum 
are 

(60) 

To get the radial eigenvectors Rnp^biP) of Eq. (001) . we 
insert the expression of the energy spectrum Eq. (go]) 
into Eq. B71) . Then, we get from Eq.® and Eq. ® : 

K,,. iy) = iy) 

(61) 

where q and p are given in Eq. ((501) . 

After inserting Eq. (ICTl) into Eq. (HOI) , we obtain 


Xuf,,Liy) 


/ 1 -l- 2 T — 2np 
\ y(l + ay) 



Xn^,Liy) 


I anp{np + |) 
\ ?/(l + ay) 


Xnp,L{y) 


(62) 


The excited state wave functions of this equation are ob¬ 
tained through Eq. (IA2I) 

xiy) = N^^_^2Ei{-np,-np - ^\-2np - + ay) 

(63) 


where Nn^^L is a normalization constant and 2 F 1 are 
hyper-geometrical functions. Therefore, according to 
the relation between hyper-geometrical functions and the 
generalized Jacobi polynomials, Eq. (4.22.1) of Ref. 
the radial wave function can be written as 


R. 


rifi,Lit) = N, 






(1 _ +p)/ 4(2 + ^)(9-tl)/4p^9/2.p/2)(^)^ 


-1 + ay 
1 +ay 


(64) 


To determine we use the usual orthogonality re¬ 

lation of Jacobi polynomials Eq. (7.391.7) of Ref. [s^. 
This leads to 

Nnii,L = (^2a'’/^+^n;3!) "" 

/ F(n ;3 + ^ + l)r( 2 n^ + ^ + l + |^) 

\r(n/3 + I + l)r(n/3 + -g^ + 2 )b(2n/9 -I- ^ + l) / 

(65) 
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a. In the case where Bp = = Brot = 1 and 

fi = 1, the wave function Eq. (1641) and the normalization 
constant Eq. (l65ll match up with Eq. (108) and Eq. 
(112) of Ref. |l4l) respectively. 

b. In the limit case a —> 0, no dependence of the 
mass on the deformation, the second-order differential 
equation Eq. (1401) must have a solution of the form 

Rn,AP) = ( 66 ) 


where b = ./ . By using this radial function in Eq. 

y -^Po 

(1401) and introducing a new variable y = , one can get 


G",,l(2/) = -' 


'l + f-26y 


(67) 


Erom Eq. m of lAM, one can define Ao(0) and so(y). 
Then, \n{y) and Sn{y) are calculated by the recurrence 
relations given in Eq. (IA4I) and the solution of this equa¬ 
tion is found through Eq. (IA2|) to be 

Gn,Ay) = {^hy) (68) 

where L denotes the Laguerre polynomials, is a 

normalization coefficient determined from the normaliza¬ 
tion condition 




F^P)dp=l 


(69) 


leading to 


Nnf,,L = 


2(26) 




np\ 


T{np + \qn^{L,k) + 1) 


1/2 


(70) 


IX. B(E2) TRANSITION RATES 


The electric quadrupole operator for axially deformed 
nuclei around 7 = 0 is given by 


rp{E2) 


tP 


COS7-b 




( 2 ) 
M,2 ' 


A 

d,- 


I sm7 


(71) 


where t is a scaling factor. The first term describes IS.K = 
0 transitions and the second is for AA" = 2 transitions. 

The B{E2) transition rates from an initial to a final 
state are given by [s^ 


B{E2-UK, 


LfKf) 


167 r 2Li + l 

( 72 ) 


and the reduced matrix element can be obtained by using 
the Wigner-Eckrat theorem [s^ 


= (L,A-,||T<g^>||L.K.) (73) 

y^2Lf + l 


The final result 


reads 


B{E2-npLn^K —^ n'pL'A^') 

= 2, K' - K\L\ L'Pn^iC.n'iC 


(74) 


with 


Inf,,L,n’p,L' = j P El ,u ^{P) Fl-P ^ dP 

= J PRL,AP)RL',n'^{P)dp (75) 

C contains the integral over 7 . For AAT = 0 cor- 

n^K,n'^K' ^ ' 

responding to transitions {g.s. ^ P ^nd 

P —>■ g.s.), the 7 —integral part reduces to the orthonor¬ 
mality condition of the 7 -wave functions : = 

S / „, ■ While for AK = 2 corresponding to transi- 

n~y,Ti^ K,K 

tions (7 —>■ g.s.,j —>• P), this integral takes the form. 

CnyK,n'^K' = J sm ^rjriyKVn'^K'l sinSjA (76) 

In the next sections, all values of B{E2) are calculated 
in units of i?(A 2 ; 2 ^ —> 0 )*'). 

X. NUMERICAL RESULTS AND DISCUSSION 

Before starting any calculations of the energy spec¬ 
tra and transition rates for the axially symmetric pro¬ 
late deformed nuclei ^^^Sm, ^^®Gd, ^’^^Yb, and 
which have been the object of Ermamatov et al. study 
[ 2 ^ and before trying to improve them within DDMF, 
we have to reevaluate the parameters of the problem 
through the corrected formulas of these nuclear char¬ 
acteristics Eqs. (j54l55D . For this purpose, we deter¬ 
mine the free parameters Bj/Bp,g, and gp from ex¬ 
perimental data of E{2A/E{2^), E{0p)/E{2^) and 
R(A 2 ; 2 + —;> 0 )''/R(A 2 ; 2 )*" —>• 0 )*') by solving a sys¬ 
tem of three nonlinear algebraic equations (Appendix B), 
while Bp/Brot is fixed to the value given in [l^. With 
the new parameters (Table [T| we have calculated the cor¬ 
rect values that Ermamatov et al. [ 2 ^ should obtain for 
the ratios A(L+j,)/A( 2 +j, ) for the ground state band, 
A(L^)/A(2+g) for /3-band and A(L+)/A(2+g ) for the 
7 -band. Here E{Lf) {i = g.s.,P,j) is the energy of the 
level characterized by the angular momentum L'^ in the 
band i and A( 2 +g ) the energy of the hrst excited level 














TABLE 1. The values of free parameters used in the calculations 


nucleus 

9 

9P 

Bp/By 

Bp ! Brot 

gi^Bp — B-y — Brot^ 

gpi^Bp — Bry — Brot) 


0.0187 

281.66 

1.36 

3.99 

0.0489 

0.357 


0.0252 

308.84 

1.53 

4.64 

0.0673 

0.884 

172Yb 

0.0064 

2469.46 

1.32 

11.14 

0.0453 

-1.909 

182W 

0.0249 

619.74 

2.01 

6.62 

0.0714 

0.512 


of the ground state band. As a qualitative test of agree¬ 
ment between the theoretical results and the experimen¬ 
tal data, we evaluated the rms differences given by 


(n-l)A( 2+)2 

where Ei (exp) is the experimental energy of the level, 
Ei(th) the corresponding theoretical value, n the max¬ 
imum number of considered levels and E{2'^) the head 
energy of the band under consideration. 

In table El we compare our results for ^^^Sm in the 
both cases Bp ^ ^ Brot (the third column with a = 

0) and Bp = B^ = B^ot (the fifth column with a = 
0) with experimental data and the data from Ref. 

One can see that our results for Bp ^ B^ ^ Brot 
agree with experimental data, particularly in /3 and 7 
bands (cr < 1 ) but are slightly different from data of 
0 . This slight discrepancy could be reduced in the 
frame of DDMF. While in the g.s. band the precision 
of our results (a > 1) is obviously affected by the energy 
value of the level L = 12 which is nearly 10% higher than 
the experimental one. From the same table, we can also 
see that the obtained values in the case Bp ^ B^ ^ Brot 
are more precise {atotai < 1) than those for which Bp = 
Bry = Brot {(^totai > !)■ For ^®®Gd (table ED) our results 
are relatively better in the 7 band for Bp 7 ^ B^ ^ Brot 
but are globally more precise than for Bp = Bj = Brot- 
Moreover, our energy spectrum for ^^^Yb given in table 
|TV]well reproduce the standard ones, particularly in the 
g.s and 7 bands with Bp ^ B-y ^ Brot unlike those of 
the case where these mass parameters are taken to be 
equal to one. On the other hand, our results for the 
nucleus (table IIVI) are more accurate (cr < 1) in 

the three bands with Bp 7 ^ B.y 7 ^ Brot than in the case 
of Bp — By — Brot • 

In order to improve the obtained numerical results, 
we recalculated the energy ratios in the framework of 
DDMF with the more elaborated formula given in Eq. 
dUl). Such a formula contains two supplementary pa¬ 
rameters, namely : a and /3o- The optimal values of both 
parameters are evaluated through rms fits of energy lev¬ 
els by making use of Eq. (1771) for each band of each 
nucleus. 

Erom tables HUTVI one can see that a fair enhancement 
of numerical results has been achieved within DDMF in 
both cases : Bp ^ By ^ Brot and Bp = By = Brot- 

Indeed, from the numerical calculations for nuclei 
^®^Sm, ^®®Gd, ^^^Yb, and we remark that the pre¬ 


cision in the case of Bp By Brot increases with the 
mass number. 

Similarly, we have also calculated transition rates 
BiE2-L+,,+2 L+J, B{E2-L+ L+,) and 

B{E2-L+ L+,,) in units of 5(^2; 2+,. ^ 0+,.) 

for the same nuclei in both cases : Bp 7 ^ By 7 ^ Brot 
and Bp = By = Brot within and out DDMF. Within 
DDMF, we have used the same optimal values of the two 
parameters a and /3o previously obtained for the energy 
ratios. 

Then in tables IVUlIXi it is clearly shown that our re¬ 
sults in the case of Bp 7 ^ By 7 ^ Brot are better than 
those with Bp = By = Brot- Indeed, in the (3 band and 
in the case of different mass coefficients with a = 0 , the 
mean difference between the theoretical value of tran¬ 
sition rate and the experimental one corresponding to 
transition 2 ^ — 0 +,, is about 1 . 8 , while in the case of 
equal mass parameters, it is equal to 20. Likewise, in the 
transition 2 ^ —4+^ , the mean difference between the 
theory and the experiment in the first case is about 15.9 
and in the second one it is equal to 152. For a 7 ^ 0, in 
the case of different mass coefficients, for the same tran¬ 
sitions 2~p —>• 0 +g and 2^ —^ Og.s.i the mean difference 
value is about 1.7 and 14, respectively, while in the case 
of equal mass coefficients it is equal to 20.3 and 152.5, 
respectively. Such a fact can also be seen in the 7 band. 
We underline here that, in the case of equal mass param- 
eters, the obtained results reproduce those of Bonatsos et 
al. [Q. The slight difference between them came from 
the fact that the Bonatsos et al. fitting calculations have 
been carried on a given number of levels which is different 
from the number we considered in our calculations. This 
is a further proof that our formulas given in Eq. (1491) and 
Eq. (IMl) respectively for the energy and the wave func¬ 
tions are more accurate than those erroneously derived 
by Ermamatov et al. [ 2 ^. Moreover, this comparison 
corroborates the fact that the mass parameter should be 
taken into account in such calculations. As it has been 
mentioned in the introduction, the Bohr Hamiltonian is 
a quite competitive method in respect to other methods 
like IBM-1 Q. To make a simple comparison between 
them, we give in tables IXlIXIl our obtained results com¬ 
pared with the available IBM-1 data. 

XI. CONCLUSION 

In this paper we have revisited all calculations per¬ 
formed in a recent work based upon inaccurate for- 
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mulas for the energy spectrum and transition rates for 
axially symmetric prolate nuclei. With the asymptotic 
iteration method we have derived the correct formulas 
for these nuclear observables. Also, we have extended our 
calculations into deformation dependent effective masses 
formalism in order to improve the numerical results. 
Moreover, we have shown the importance of the mass 
parameter to be introduced in numerical calculations un¬ 
like what it has been done by other authors who have 
neglected the important role played by this parameter in 
such calculations. Through a comparison with IBM-1, 
the Bohr Hamiltonian with mass parameters has proved 
to be more accurate. 


Appendix A: Asymptotic Iteration Method (AIM) 


The asymptotic iteration method is proposed to 
solve the second-order homogeneous differential equation 
of the form 


y''{x) = Xo{x)y'{x) -f so{x)y{x) (Al) 

where the variables Aq and sq are sufficiently differen¬ 
tiable. 

The differential equation m has a general solution (25j | 


y{x)=exp(^ — J a{xi)dxi^ C 2 

-\-Ci J exp ^ y [Ao(a; 2 ) + 2 Q;(a; 2 )]da; 2 ^c!a;i (A2) 


the energy eigenvalues are then computed by means of 
the following termination condition [25| 

S — SnXn—l XnSji—l — 0 (-^^) 

Appendix B: Formulas used for the calculations of 

the B{E2) 


In this appendix we present the expressions used for 
calculations of the transition probabilities B{E2) : 

H(A2;T;+. ^£+,.) 

H(A2;2+,. ^0+..) 

= 5(CfPo2o)^ 

^ / r[0.5(go(T',0)+go(£,0)) + 1.5] y 
"" r[0.5(go(2,0) + go(0,0))+1.5] ) 

^ r[go(2,0) + l]r[go(0,0) + l] 

r[go(T',0) + l]r[go(T,0) + l] 


B{E2;L'+^L+J 
B{E2;2+s. ^ 0 +«.) 

_ ^/r'LO \2 

— ^Wl'0201 

^ / r[0.5(go(L',0)+go(L,0)) +1.5] \ 
"" r[0.5(go(2,0) + go(0,0))+1.5] ) 

^ r[go(2,0) + l]r[go(0,0) + l] 

"" r[go(T',0) + l]r[go(T,0) + l] 

^ (go(L',0)-go(L,0)-l)2 

QoiL', 0) -I- 1 


(B2) 


If we have n > 1, for sufficiently large n, a{x) values 
can be obtained 


Sn{x) _ Sn-l{x) 
Xn{x) Xn-l{x) 


(A3) 


with the sequences 


A„(x) =A'„_i(a;) -I- s„_i(x) -I- Ao(x)A„_i(a;) (A4a) 

s„(a;) =s'^_j^{x) + so{x)Xn-iix), n = 1,2,3... 

(A4b) 


BiE2;L'+^L+J 
B{E2-,2+s. ^0+«,) 

= 5g(CfPg2o)^ 

^ / r[0.5(go(L',0)+go(L,0)) +1.5] \ 
"" r[0.5(go(2,0) + go(0,0))+1.5] ) 

^ r[go(2,0) + l]r[go(0,0) + l] 

"" r[go(T',0) + l]r[go(T,0) + l] 

where C^!l) 2 o Clebsch-Gordan coefficients. 


(B3) 
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TABLE 11. The comparison of the theoretical predictions of energy levels Eq. (|49l) of the ground state band, the P and 7 bands 
normalized to the energy of the first excited state using the parameters given in Table U for for this work with 

those from Ref. [T^ and experimental values taken from Ref. [^. Po and a indicate the position of the minimum of Davidson 
potential Eq. (1421) and the deformation dependence of the mass Eq. (1431) respectively, while a is the quality measure Eq. (1771) . 


L 

exp 

Bp ^ B^l ^ Brot 

Bp = 

S'j — Brot 

Ref. [IJ 

a = 0 

DDM 

a — 0 

DDM 

g.s. 







4 

3.26 

3.31 

3.31 

3.25 

3.27 

3.28 

6 

6.63 

6.89 

6.89 

6.59 

6.68 

6.76 

8 

11.01 

11.65 

11.65 

10.82 

11.07 

11.28 

10 

16.26 

17.52 

17.52 

15.75 

16.29 

16.65 

12 

22.27 

24.41 

24.41 

21.22 

22.24 

22.68 

CT 


1.289 

1.289 

0.592 

0.043 

0.320 

a 



0.0000 


0.0483 


Po 



22.41 


0.54 


Pi 







0 

13.40 

13.40 

13.03 

13.40 

13.13 

13.40 

2 

14.37 

14.40 

14.03 

14.40 

14.14 

14.40 

4 

16.32 

16.71 

16.35 

16.65 

16.41 

16.68 

6 

19.23 

20.29 

19.94 

19.99 

19.78 

20.16 

8 


25.05 

24.73 

24.22 

24.08 

24.68 

10 


30.93 

30.64 

29.15 

29.11 

30.05 

12 


37.81 

37.58 

34.62 

34.73 

36.08 

a 


0.651 

0.501 

0.479 

0.384 

0.576 

a 



0.0335 


0.0039 


Po 



0.88 


0.88 


71 







2 

17.56 

17.56 

18.01 

17.56 

18.67 

17.56 

3 

18.77 

18.47 

18.97 

18.28 

19.49 

16.56 

4 

20.30 

19.68 

20.26 

19.23 

20.57 

19.87 

5 

22.01 

21.18 

21.86 

20.39 

21.91 

21.48 

6 

23.73 

22.96 

23.78 

21.77 

23.50 

23.38 

7 

26.27 

25.03 

26.01 

23.33 

25.33 

25.53 

8 


27.36 

28.56 

25.08 

27.40 

27.93 

9 


29.95 

31.41 

26.99 

29.70 

30.55 

10 


32.79 

34.58 

29.05 

32.21 

33.36 

11 


35.88 

38.06 

31.26 

34.94 

36.33 

12 


39.19 

41.84 

33.60 

37.87 

39.44 

13 


42.73 

45.93 

36.06 

40.99 

42.65 

CT 


0.812 

0.260 

1.817 

0.743 

1.097 

a 



0.0199 


0.0086 


Po 



11.53 


1.67 




0.895 

0.874 

1.153 

0.754 

0.728 

a 



0.0219 


0.0054 


Po 



1.07 


1.60 
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TABLE III. The comparison of the theoretical predictions of energy levels Eq. (1491) of the ground state band, the P and 7 bands 
normalized to the energy of the first excited state i5(2^s,) using the parameters given in Table|T]for ^®®Gd for this work with 
those from Ref. [3 s^nd experimental values taken from Ref. [^. /3o and a indicate the position of the minimum of Davidson 
potential Eq. (1421) and the deformation dependence of the mass Eq. (1431) respectively, while cr is the quality measure Eq. (1771) . 




Bp ^ B^l ^ Brot 

Bp = 

.S'j — Brot 


L 

exp 

a = 0 

DDM 

a — 0 

DDM 

Ref. [ij 

g.s. 

4 

3.24 

3.31 

3.30 

3.23 

3.25 

3.29 

6 

6.57 

6.87 

6.82 

6.49 

6.60 

6.76 

8 

10.84 

11.62 

11.50 

10.55 

10.87 

11.22 

10 

15.91 

17.44 

17.30 

15.22 

15.92 

16.51 

12 

21.63 

24.25 

24.18 

20.34 

21.62 

22.38 

cr 


1.575 

1.492 

0.747 

0.025 

0.470 

a 



0.0600 


0.0217 


Po 



60.00 


1.11 


Pi 

0 

11.79 

11.79 

9.93 

11.79 

11.79 

11.79 

2 

12.69 

12.79 

10.94 

12.79 

12.79 

12.79 

4 

14.68 

15.10 

13.29 

15.02 

15.02 

15.08 

6 

17.30 

18.66 

16.93 

18.28 

18.28 

18.54 

8 

20.76 

23.41 

21.82 

22.34 

22.34 

23.01 

10 

24.94 

29.23 

27.89 

27.01 

27.01 

28.30 

12 

30.43 

36.04 

35.09 

32.13 

32.13 

34.17 

a 


3.135 

2.585 

1.339 

1.339 

2.311 

a 



0.0230 


0.0000 


Po 



2.77 


43.81 


71 

2 

12.97 

12.97 

12.97 

12.97 

13.75 

12.97 

3 

14.02 

13.90 

13.90 

13.70 

14.58 

13.96 

4 

15.22 

15.12 

15.12 

14.66 

15.67 

15.27 

5 

16.93 

16.64 

16.64 

15.83 

17.02 

16.88 

6 

18.47 

18.45 

18.45 

17.20 

18.61 

18.76 

7 

20.79 

20.54 

20.54 

18.75 

20.44 

20.90 

8 

22.60 

22.90 

22.90 

20.47 

22.50 

23.27 

9 

25.28 

25.51 

25.51 

22.34 

24.78 

25.85 

10 

27.44 

28.38 

28.38 

24.35 

27.28 

28.60 

11 

30.19 

31.48 

31.48 

26.49 

29.97 

31.50 

12 

32.84 

34.80 

34.80 

28.74 

32.86 

34.52 

13 

35.67 

38.34 

38.34 

31.08 

35.94 

37.63 

cr 


1.122 

1.122 

2.725 

0.390 

0.982 

a 



0.0000 


0.0219 


Po 



32.35 


1.42 




1.897 

1.866 

2.029 

1.008 

1.379 

a 



0.0499 


0.0582 


Po 



0.98 


0.80 
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TABLE IV. The comparison of the theoretical predictions of energy levels Eq. (gSll of the ground state band, the P and 7 bands 
normalized to the energy of the first excited state using the parameters given in TableUfor ^^^Yb for this work with 

those from Ref. [3 s^nd experimental values taken from Ref. [^. /3o and a indicate the position of the minimum of Davidson 
potential Eq. (1421) and the deformation dependence of the mass Eq. (1431) respectively, while cr is the quality measure Eq. (17711 . 




Bp ^ B^l ^ Brot 

Bp = 

.S'j — Brot 


L 

exp 

a = 0 

DDM 

a — 0 

DDM 

Ref. [IJ 

g.s. 

4 

3.29 

3.33 

3.30 

3.25 

3.28 

3.32 

6 

6.84 

6.96 

6.84 

6.58 

6.78 

6.91 

8 

11.54 

11.87 

11.54 

10.79 

11.45 

11.71 

10 

17.34 

18.02 

17.33 

15.69 

17.26 

17.65 

12 

24.14 

25.36 

24.15 

21.12 

24.21 

24.64 

cr 


0.719 

0.087 

1.766 

0.076 

0.309 

a 



0.0036 


0.0800 


Po 



40.09 


1.98 


Pi 

0 

13.20 

13.20 

11.68 

13.20 

13.20 

13.20 

2 

14.15 

14.20 

12.69 

14.20 

14.20 

14.20 

4 

16.34 

16.53 

15.02 

16.45 

16.45 

16.52 

6 

19.53 

20.16 

18.67 

19.78 

19.78 

20.11 

8 

23.54 

25.07 

23.62 

23.99 

23.99 

24.93 

10 

28.10 

31.22 

29.83 

28.89 

28.89 

30.89 

12 

33.11 

38.56 

37.27 

34.32 

34.32 

37.91 

a 


4.593 

2.135 

0.628 

0.628 

2.311 

a 



0.0418 


0.0000 


Po 



1.74 


1.60 


71 

2 

18.63 

18.63 

18.71 

18.63 

19.08 

18.63 

3 

19.68 

19.59 

19.69 

19.33 

19.89 

19.63 

4 

21.06 

20.87 

20.98 

20.26 

20.95 

20.95 

5 

22.60 

22.47 

22.60 

21.39 

22.28 

22.60 

6 


24.38 

24.54 

22.73 

23.87 

24.56 

7 


26.60 

26.80 

24.26 

25.73 

26.83 

8 


29.12 

29.37 

25.96 

27.84 

29.39 

9 


31.95 

32.26 

27.83 

30.20 

32.24 

10 


35.07 

35.46 

29.85 

32.83 

35.37 

11 


38.48 

38.97 

32.01 

35.70 

38.76 

12 


42.17 

42.78 

34.29 

38.83 

42.20 

13 


46.15 

47.91 

36.70 

42.22 

46.28 

cr 


0.121 

0.065 

0.862 

0.347 

0.070 

a 



0.0075 


0.0600 


Po 



17.10 


2.19 




1.719 

1.413 

1.067 

3.458 

1.495 

a 



0.0010 


0.0100 


Po 



11.12 


90.01 
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TABLE V. The comparison of the theoretical predictions of energy levels Eq. (|49l) of the ground state band, the P and 7 bands 
normalized to the energy of the first excited state lf(2g,s,) using the parameters given in Table U for for this work with 

those from Ref. [3 s^nd experimental values taken from Ref. [^. /Jq and a indicate the position of the minimum of Davidson 
potential Eq. (1421) and the deformation dependence of the mass Eq. (1431) respectively, while cr is the quality measure Eq. (17711 . 


L 

exp 

Bp ^ B^l ^ Brot 

Bp = 

.S'j — Brot 

Ref. [IJ 

a = 0 

DDM 

a — 0 

DDM 

g.s. 







4 

3.29 

3.32 

3.29 

3.22 

3.29 

3.30 

6 

6.80 

6.91 

6.78 

6.45 

6.78 

6.81 

8 

11.44 

11.71 

11.40 

10.47 

11.40 

11.41 

10 

17.12 

17.65 

17.07 

15.05 

17.07 

16.94 

12 

23.72 

24.64 

23.77 

20.07 

23.76 

23.2 

cr 


0.548 

0.042 

2.161 

0.037 

0.276 

a 



0.0335 


0.0470 


Po 



18.49 


1.25 


Pi 







0 

11.36 

11.36 

11.36 

11.36 

11.36 

11.36 

2 

12.57 

12.36 

12.36 

12.36 

12.36 

12.36 

4 

15.10 

14.68 

14.68 

14.58 

14.58 

14.66 

6 


18.26 

18.26 

17.81 

17.81 

18.17 

8 


23.07 

23.07 

21.83 

21.83 

22.77 

10 


29.01 

29.01 

26.41 

26.41 

28.30 

12 


36.00 

36.00 

31.43 

31.43 

34.57 

a 


0.335 

0.335 

0.395 

0.395 

0.345 

a 



0.0000 


0.0000 


Po 



53.35 


52.55 


71 







2 

12.21 

12.21 

12.42 

12.21 

12.76 

12.21 

3 

13.31 

13.16 

13.19 

12.94 

13.55 

13.21 

4 

14.43 

14.41 

14.46 

13.89 

14.60 

14.52 

5 

16.24 

15.97 

16.03 

15.05 

15.90 

16.14 

6 

17.70 

17.83 

17.90 

16.41 

17.42 

18.05 

7 


19.71 

20.07 

17.95 

19.17 

20.24 

8 

22.61 

22.41 

22.52 

19.65 

21.13 

22.68 

9 


25.11 

25.26 

21.50 

23.28 

25.35 

10 


28.08 

28.26 

23.48 

25.62 

28.24 

11 


31.31 

31.53 

25.58 

28.14 

31.31 

12 


34.78 

35.06 

27.79 

30.83 

34.55 

13 


38.49 

38.83 

30.10 

33.68 

37.92 

cr 


0.168 

0.158 

0.935 

0.380 

0.194 

a 



0.0538 


0.0215 


Po 



0.95 


1.07 




0.358 

0.357 

1.369 

1.019 

0.240 

a 



0.0000 


0.0054 


Po 



50.35 


2.16 
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TABLE VI. The comparison of the theoretical predictions of B{E2) in units of 5(172; 2^^ —y O^s.) using the parameters 
given in Table H] for ^®"^Sm in this work with those from Ref. [I^ and experimental values. 




exp 

Bp ^ Bry ^ B-rot 

a = 0 DDM 

Bp = 

a = 0 

B'y - Brot 

DDM 

Ref. [18] 

B{E2-L^,s 

.+2->--£'q.a.) 







B{E2-,2^, 

.s. ^O^.s. ) 







4.S. 

2+ 

^g.s. 

1.40(5) 

1.44 

1.44 

1.47 

1.47 

1.44 

Ks. 

^ts. 

1.67(7) 

1.61 

1.61 

1.70 

1.70 

1.61 

8^., ^ 

^Is. 

1.83(11) 

1.72 

1.72 

1.90 

1.90 

1.72 

10^.. — 

> 8+., 

1.81(11) 

1.81 

1.81 

2.10 

2.10 

1.82 

12+., — 

> lot. 


1.90 

1.90 

2.30 

2.30 

1.91 


B(£/2;Z/^- yLg.s .) 

fl(B2;2+«.->0+.s.) 


X 10® 


2+ 

i 

"I 

6« 


Oj-s 

^t.s 

4+ 

^g.s 

2 "*" 

^g.s 

4+ 

^Q.S 

Gts 

2 + 

4+ 

^Q.S 

Ks 

8^. 


B(E2;L^ >Lg,s.) ^ 

B(E2-2+^_^0+s.) 


2 + 

2j 

2 + 

■^7 

q+ 


^7 

4"^ 

^7 

4'^ 

^7 


5 

5 


+ 

7 

+ 

X. 


Og.s 

2^. 

4+ 

2+ 

•^g.s 

4+ 

^g.s 

2+ 

4^ 

^q.s 

^t-s 

4+ 

^q.s 


6 


+ 

.g.s 


5.4(13) 


25(6) 


18.4(29) 

3.9(6) 


6.5 

6.4 

23.3 

24.3 

6.7 

5.5 

5.5 

20.0 

21.3 

5.6 

3.1 

3.1 

12.4 

13.7 

2.9 

12.9 

11.9 

46.4 

47.9 

13.3 

11.7 

11.6 

42.3 

43.6 

12.1 

11.5 

11.4 

41.5 

42.9 

11.9 

59.7 

59.0 

216.5 

221.5 

61.7 

42.2 

35.6 

152.0 

154.4 

43.8 

48.5 

47.7 

169.9 

171.8 

51.3 

57.4 

56.2 

191.6 

193.1 

62.8 


18.4 

14.8 

46.5 

48.7 

18.4 

26.5 

21.3 

68.6 

71.4 

26.2 

1.4 

1.1 

3.7 

3.8 

1.3 

33.1 

29.3 

85.0 

88.8 

32.8 

13.5 

11.8 

36.5 

37.6 

13.0 

11.0 

9.7 

28.0 

29.4 

11.0 

33.0 

28.7 

88.6 

91.7 

32.1 

2.9 

2.4 

8.4 

7.8 

2.7 

29.8 

25.8 

79.4 

82.65 

29.3 

17.6 

14.8 

49.9 

51.0 

16.5 
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TABLE VII. The comparison of the theoretical predictions of B{E2) in units of B{E2-,2'^^ O^g.) using the parameters 

given in Table H] for in this work with those from Ref. [I^ and experimental values. 


Bf^ ^ B'y ^ Brot Bf^ — B^ — Brot 


exp 

a = 0 

DDM 

a = 0 

DDM 

Ref. [1 

1.41(5) 

1.44 

1.37 

1.48 

1.48 

1.44 

1.58(6) 

1.61 

1.42 

1.73 

1.73 

1.61 

1.71(10) 

1.72 

1.38 

1.96 

1.95 

1.73 

1.68(9) 

1.82 

1.32 

2.19 

2.18 

1.83 

1.60(16) 

1.91 

1.26 

2.43 

2.42 

1.93 

3.4(3) 

6.1 

4.4 

24.7 

24.7 

6.3 


4.6 

3.9 

19.5 

19.5 

4.7 


2.1 

1.9 

11.0 

11.0 

0.9 

18(2) 

12.6 

4.2 

51.9 

51.9 

13.0 


11.5 

9.4 

47.3 

47.3 

11.8 


11.2 

3.6 

46.4 

46.4 

11.6 


60.5 

48.9 

251.1 

251.1 

62.5 

22(2) 

44.2 

1.9 

181.4 

181.4 

46.0 


52.0 

41.2 

204.9 

204.9 

53.3 


62.4 

48.9 

230.6 

230.6 

68.9 


B ( g 2 ; I .„,.,+2 - 


S(B2;2j.s 


L 

Oj.s.) 


d”*" - 

^q.s. 
6^s. - 
8^- - 
lOs+s. 

12j-. 


■ 2 "*" 
^g.s. 

■ 

^q.s. 

6^- 

8+,. 

10^.. 


B(E2-,L0- 

B(E2-,2+^, 


X 10= 


2+ 

i 

6| 

2| 

6| 

0| 

6t 


0^. 

2 "*' 

^9-- 

■ 4"^ 

9-- 

■ 2'^ 


6+., 

■2+ 

■ 

^q.t 

6+, 


B(E2-,L0 yLg.s.) ^ 
S(B2;2+«.->0+3.) 


2 + 

2 + 

2 + 

3+ 

3+ 

4"^ 

4"^ 

4"^ 

^7 

K+ 



Og... 

25.0(8) 

25.0 

25.0 

65.9 

68.3 

25.0 


38.7(13) 

36.1 

36.1 

97.7 

100.6 

35.5 

4+ 

^9.s. 

4.1(2) 

1.8 

1.8 

5.3 

5.3 

1.8 

2+ 

^9.s. 

39.0(75) 

44.9 

44.9 

121.1 

125.2 

44.6 


27.2(35) 

18.3 

18.3 

52.3 

53.3 

17.7 

2+ 

^9.s. 

9.6(27) 

14.9 

14.9 

39.9 

41.5 

14.9 

4+ 

9'S. 

53.16(16) 

44.9 

44.9 

127.2 

130.3 

43.6 

6g.s. 


4.0 

4.0 

12.1 

11.1 

3.7 

4.S. 

43(43) 

40.5 

40.5 

114.1 

117.6 

39.8 

Gt.s. 

59(59) 

23.9 

23.9 

72.1 

72.8 

22.4 
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TABLE VIII. The comparison of the theoretical predictions of B{E2) in units of B{E2-2'^^ —> O^ ,, ) using the parameters 
giveir in Table H] for in this work with those from Ref. [l^ and experimental values. 




exp 

Bp ^ Bry ^ B-rot 

a = 0 DDM 

Bp = 

a = 0 

B'y - Brot 

DDM 

Ref. [18] 

B{E2-L^,s 

.+2->--£'q.a.) 







B{E2-,2^, 

.s. ^O^.s. ) 







4.S. 

2+ 

^g.s. 

1.42(10) 

1.43 

1.37 

1.47 

1.34 

1.43 

Ks. 

^ts. 

1.51(7) 

1.59 

1.41 

1.70 

1.36 

1.59 

8^., ^ 

^Is. 

1.89(19) 

1.67 

1.36 

1.90 

1.31 

1.67 

10^.. — 

> 8+., 

1.77(11) 

1.74 

1.29 

2.11 

1.26 

1.74 

12+., — 

> lot. 


1.79 

1.31 

2.32 

1.22 

1.79 


B(£/2;Z/^- yLg.s .) 

fl(B2;2+«.->0+.s.) 


X 10® 


2+ 

i 

"I 

6« 


Oj-s 
■ 4+ 

^g.s 

. 2~^ 
^g.s 

4+ 

^Q.S 

Gts 

2+ 

4+ 

^Q.S 

Ks 

8^. 


B(E2-,L~y- 


S(£;2;2^. 


^g-s-) 

.n+ 'i 


2"^ 

2 + 

2 + 

3+ 

3+ 

4"^ 

4"^ 

^7 

4"^ 

^7 

5+ 

5+ 


0 , 


+ 

g.s 

2+ 

4+ 

2+ 

•^g.s 

4+ 

^g.s 

2+ 

4"^ 

^Q.S 

6g.s 

4+ 

^q.s 

6g.. 


1.1(1) 

2.4 

3.6 

23.5 

23.5 

2.4 


2.0 

2.1 

19.9 

19.9 

2.0 


1.1 

1.2 

12.3 

12.3 

1.0 


4.8 

4.6 

47.1 

47.1 

4.17 


4.4 

4.3 

42.8 

42.8 

4.4 


4.3 

4.2 

42.1 

42.1 

4.3 


22.3 

21.4 

220.3 

220.3 

22.5 

12(1) 

15.9 

6.3 

155.1 

155.1 

16.0 


18.4 

17.1 

173.7 

173.7 

18.8 


22.1 

20.3 

195.8 

195.8 

22.8 

X 10® 






6.3(5) 

6.3 

3.7 

42.3 

23.4 

6.3 


9.0 

5.5 

62.6 

34.9 

9.0 

0.60(5) 

0.5 

0.3 

3.4 

1.8 

0.4 


11.3 

10.3 

77.5 

43.0 

11.2 


4.6 

4.2 

33.4 

18.1 

4.4 

33(24) 

3.8 

3.4 

25.5 

14.1 

3.8 


11.1 

10.1 

81.1 

43.8 

11.0 


1.0 

0.9 

7.7 

3.6 

0.9 


10.1 

9.1 

72.5 

39.0 

10.0 


5.8 

5.2 

45.9 

23.0 

5.7 
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TABLE IX. The comparison of the theoretical predictions of B{E2) in units of B{E2\2'^^ —> O^.s.) using the parameters 
given in Table H] for in this work with those from Ref. [I^ and experimental values. 




exp 

Bp ^ Bj ^ B-rot 

a = 0 DDM 

Bp = 

a = 0 

B'y - Bfot 

DDM 

Ref. [18] 

B{E2-L^,s 

.+2->--£'q.a.) 







B(B2;2j 









2+ 

^g.s. 

1.43(8) 

1.44 

1.36 

1.49 

1.48 

1.44 

6^s. ^ 

^ts. 

1.46(6) 

1.60 

1.40 

1.74 

1.72 

1.60 

8^.. ^ 

^Is. 

1.53(10) 

1.71 

1.35 

1.98 

1.92 

1.71 

10^.. — 

> 8+., 

1.48(9) 

1.79 

1.28 

2.22 

2.11 

1.80 

12+., — 

> lot. 


1.87 

1.23 

2.48 

2.30 

1.88 


- yLg.s .) 


X 10® 


2 + 

i 

"I 

6« 


Oj.s 

^ts 


-^g.s 

Gj.s 

2ts 

4+ 

^Q.S 

Ks 

»t.s 


B{E2-,L^- 


B(£;2;2^. 


^g-s-) 

.n+ ^ 


2"^ 

2 + 

2 + 

3+ 

3+ 

4'^ 

4"^ 

^7 

4"^ 

^7 

5+ 

5+ 


0 , 


+ 

g.s 

0 + 
^g.s 

4+ 

2+ 

•^g.s 

4+ 

2+ 

^g.s 

4+ 

^Q.S 

^t-s 

4+ 

^q.s 

69 .. 


6.6(10) 

4.4 

4.4 

25.1 

25.1 

5.2 


3.2 

3.2 

19.2 

19.2 

2.1 


1.4 

1.4 

10.5 

10.5 

0.1 

4.6(6) 

9.3 

9.3 

53.6 

53.6 

13.7 


8.4 

8.4 

48.8 

48.8 

12.5 


8.3 

8.3 

47.9 

47.9 

12.2 


44.9 

44.9 

262.3 

262.3 

77.0 

13(1) 

33.1 

33.1 

191.1 

191.1 

64.1 


39.4 

39.4 

216.3 

216.3 

81.3 


47.81 

47.81 

243.2 

243.2 

101.6 

X 10® 






24.8(8) 

24.8 

24.8 

70.2 

72.1 

24.8 

49.2(13) 

35.7 

35.7 

104.2 

106.5 

35.3 

0.2(2) 

1.8 

1.8 

5.6 

5.7 

1.8 


44.5 

44.6 

129.2 

132.5 

44.2 


18.1 

18.1 

56.0 

56.8 

17.6 

17.2(17) 

14.8 

14.8 

42.6 

43.8 

14.8 

75.9(73) 

44.3 

44.3 

136.3 

138.7 

43.3 


3.9 

3.8 

13.0 

11.8 

3.7 


40.0 

40.1 

122.2 

125.0 

39.5 


23.4 

23.4 

77.5 

78.1 

22.3 
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TABLE X. The comparison of the theoretical predictions of E{Lf) [i = g.s.,7-band) normalized to the energy of the first 
excited state using the parameters given in TableUfor ^^"‘Sm and in this work with those from IBM-1 Ref. 

and experimental values Ref. [dOll . 




164 

Sm 



18^ 

W 


L 

exp 

a = 0 

DDM 

IBM-1[^ 

exp 

a = 0 

DDM 

IBM-l[d2 

g.s. 









d 

3.26 

3.31 

3.31 

3.19 

3.29 

3.32 

3.29 

3.33 

6 

6.63 

6.89 

6.89 

7.33 

6.80 

6.91 

6.78 

6.95 

8 

11.01 

11.65 

11.65 

12.dd 

11.dd 

11.71 

11.dO 

12.00 

10 





17.12 

17.65 

17.07 

18.33 

(7 


0.d90 

0.d90 

1.127 


0.350 

0.039 

0.775 

Pi 









0 

13.d0 

13.d0 

13.03 

Id.Od 

11.36 

11.36 

11.36 

ll.dl 

2 

ld.37 

Id.dO 

ld.03 

ld.78 

12.57 

12.36 

12.36 

ll.d6 

d 

16.32 

16.71 

16.35 

17.31 

15.10 

ld.68 

ld.68 

13.81 

6 

19.23 

20.29 

19.9d 

17.66 


18.26 

18.26 

17.50 

(J 


0.651 

0.501 

1.158 


0.332 

0.332 

1.20d 

71 









2 

17.56 

17.56 

18.01 

18.53 

12.21 

12.21 

12.2d 

12.dl 

3 

18.77 

18.d7 

18.97 

18.97 

13.31 

13.16 

13.19 

12.d7 

d 

20.30 

19.68 

20.26 

21.72 

ld.d3 

Id.dl 

ld.d6 

ld.9d 

5 

22.01 

21.18 

21.86 

2d.l2 

16.2d 

15.97 

16.03 

15.d8 

6 

23.73 

22.96 

23.78 


17.70 

17.83 

17.90 

18.69 

(J 


0.623 

0.298 

1.576 


0.168 

0.158 

0.801 


TABLE XL The comparison of the theoretical predictions of B(E2) in units of B{E2\ 2t„ —> O^s.) using the parameters 
given in TableUfor and in this work with those from IBM-1 Ref. (dll. and experimental values Ref. [d^ . 




154 

Sm 



182 

W 



exp 

a = 0 

DDM IBM-1 

exp 

a = 0 

DDM IBM-1 

S(£j2,Ijp.s. +2 -) 

S(B2;2+«,-^0+,.) 









4^.s. —t 2+,. 

6+ — 4 + 

'-'g.s. ^ ^9-s- 

l.d0(5) 

1.67(7) 

l.dd 

1.61 

l.dd 

1.61 

1.35 

1.53 

l.d3(8) 

l.dd 

1.36 

1.33 

B(B2;2+«,->0+«,) 









2+ _ \ O’*" 

^ ^g s. 

5.d(13) 

6.5 

6.d 

7.78 

6.6(10) 

d.d 

d.d 

113.0 

dt — ^ 2+,. 


5.5 

5.5 

29.57 


3.2 

3.2 

56.3 

2+ _ \ 2'’" 

^ ^g.s. 


12.9 

12.8 

15.33 

d.6(6) 

9.3 

9.3 

d8.d 

2+ _ y 4+ 

^ %.s. 

25(6) 

d2.2 

dl.d 

30.67 

13(1) 

33.1 

33.1 

3.7 


B{E2-,Lf) - tLg.s.) 

B(B2;2+«,->0+«,) 


X 10® 


2+ _ 

-^0+,. 

18.4(29) 

18.4 

16.5 

16.43 


2+ _ 
^7 

-^2+,. 


26.5 

23.6 

2d.l0 


2+ — 
Z7 

4+,. 

3.9(6) 

l.d 

1.2 

0.99 

0.2(2) 1.8 1.8 153.6 


























